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Abstract 

We have computed the SU{2) Low Energy Constant ^5 and the mass splitting between charged and 
neutral pious from a lattice QCD simulation of = 2 + 1 flavors of Domain Wall Fermions at a scale of 
a^^ = 2.33 GeV. Relating /s to the S parameter in QCD we obtain a value of S{mH — 120 GeV) = 0.42(7), 
in agreement with previous determinations. Our result can be compared with the value of S from electroweak 
precision data which constrains strongly interacting models of new physics like Technicolor. This work in 
QCD serves as a test for the methods to compute the S parameter with Domain Wall Fermions in theories 
beyond the Standard Model. We also infer a value for the pion mass splitting in agreement with experiment. 

PACS numbers: 11. 15. Ha, ll.30.Rd, 12.38.Gc 12.39.Fe 



'Electronic address: |j wenneke@ph.ed.ac.uk 



I. INTRODUCTION 



Models of dynamical electroweak symmetry breaking (DEWSB) like Technicolor are possible 
extensions of the Standard Model, and many proposals have been put forward for a strongly- 
interacting sector beyond the Standard Model since the original proposals in Refs. [U |2]. Viable 
models of DEWSB must satisfy the constraints that follow from electroweak precision data [3l S]. 
These constraints put severe limitations on Technicolor candidates, and QCD-like theories naively 
rescaled to the electroweak scale are already ruled out, see e.g. Refs. [5, 6j for recent reviews. The 
constraints are nicely encoded in bounds for the value of the S parameter introduced in Refs. [31 [7]. 

Theories near an infrared fixed point (IRFP) have been advocated as promising candidates 
for DEWSB. Preliminary attempts at estimating analytically the S parameter for these theories 
suggest that it is much reduced compared to QCD-like theories [H [9]. Unfortunately these com- 
putations are not based on first principles, and have to rely on assumptions that are difficult to 
control. Early studies focused on finding evidence for an IRFP in theories with a large number 
of flavors in the fundamental representation of the gauge group following the seminal example in 
Ref. [lOj. More recently, novel models have been proposed based on smaller numbers of fermions in 
higher-dimensional representations [11]. Several phenomenological scenarios have been proposed 
which build upon these ideas [12l [13] . 

Lattice simulations are now in a position to address these questions from first principles, so that 
the difficulties in dealing with the non-perturbative dynamics can be dealt with in a systematic 
way. The existence of IRFPs has been investigated in theories with fundamental fermions |14| [TS] 
[ISl [El EHl [13 ED] and higher representations [2ll|22l[23llSl|25l[Sll2Zl[Mll29ll3Ql[lI]. These 
preliminary studies have mapped out the space of bare lattice parameters and are starting to study 
the spectrum of the candidate theories, and the flow of renormalized couplings. First results hint 
towards an interesting landscape of theories that could exhibit scale invariance at large distances. 

Computing the S parameter in these theories from first principles is an important ingredi- 
ent in trying to build successful phenomenological models. The S parameter is obtained on the 
lattice from the form factors that appear in the momentum-space VV-AA correlator, as defined 
below in Section |IT} Chiral symmetry plays an important role in guaranteeing the cancellation 
of power-divergent singularities when computing the above correlator. Hence lattice formulations 
that preserve chiral symmetry at finite lattice spacing are particularly well-suited for these studies. 
QCD is the ideal testing ground to develop and test the necessary lattice technology. In this work, 
we compute the S parameter in QCD with nj = 2 + 1 fiavors of Domain Wall fermions (DWF). 
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Our study closely follows the procedure described in Ref. [32] where the S parameter (or equiva- 
lently the 5'f7(3) Low Energy Constant Lio) was first computed from vacuum polarisation functions 
(VPFs) using overlap fermions. We adopt the method, apply it to Domain Wall Fermions, and 
widen its scope by using conserved currents and larger physical volumes. 

Section In] contains a short account of computational methods for the VPF on the lattice. 



Section III addresses the topics of power divergences and residual chiral symmetry breaking. In 
Section IV we present the numerical results obtained on the gauge configurations produced by the 
RBC and UKQCD collaborations using Domain Wall Fermions. Section [V] contains our results for 
the pion mass splitting. A discussion of the numerical results and a short conclusion can be found 
in Section IVTl 

II. VACUUM POLARISATION FUNCTIONS 

Domain Wall Fermions are a five dimensional formulation of lattice QCD with an approximate 
chiral symmetry |33| IM] [35] . The residual explicit breaking of chiral symmetry appears in the Ward 
identities of the theory as terms proportional to the so-called residual mass m^es- The conserved 
vector and axial currents form a multiplet under this approximate lattice chiral symmetry. 

The basic observables in this work are vacuum polarisation functions of the vector and the axial 
vector current. They are defined as current-current two-point functions in momentum space, 

n;:.(g) = 5]e'^-(0|V^(x)K(0)|0), (1) 

X 

ii^M^T.^''"'(^\M^)Mom, (2) 

X 

where and are the conserved vector and axial currents and and are the corresponding 
local currents. We consider local-conserved correlators from a new set of point source propagators 
with up to two units of spatial momentum. The definition of the conserved currents can be 
found in Ref. [35]. Since the Fourier transform in Eqs. (lj2) includes x = 0, power-divergent 



contributions can arise. Due to lattice chiral symmetry these divergences cancel in the difference of 
the vector and axial vector correlators, if conserved currents are used. Note that similar observables 
are used to compute hadronic contributions to the anomalous magnetic moment of the muon in 
Refs. |36| 157] l38] . where one can find a more detailed discussion of the renormalization of the 
correlators. A similar cancellation was pointed out in Ref. [32], where overlap fermions have been 
used. 
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Following Ref. 

transverse part, 



we decompose the difference 11 



V-A 



H'^j^ into a longitudinal and a 



n 



V-A 



(3) 



For each momentum we average all components of Ilfj^u which contribute to only one of the "polar- 
isations" . For example for momenta with one non-zero direction we have g2n(0) = n^"^ and 
2]-[(i) = 1 ^ n^-^ 

We use SU{2) Chiral Perturbation Theory (ChPT) to fit to our data to extract the Low Energy 
Constant Z^. Since we are using 2-1-1 flavor lattices the low-energy constants will implicitly depend 
on the strange quark mass. The more common notation in the literature is to use the SU{3) LEG 
L'lQ which is related to by 



rr jr 

^10 — '5 



1 



log 



m 



K 



3847r2 V /i 
= 11-2,- 10"^ for II = m 

The ChPT result for H^^) can be found in Ref. [39]: 

1 1\ 1 



+ 1 



(4) 



p- 



g^n^^) (m. 
with J{a 



1 



487r2 



, a - 1 
o-log — — r + 2 



3 

and a 



(5) 
(6) 



167r2 V" a + 1 

With m^, known from fits to pseudoscalar and Aq correlators (Table |l]) there is only one free 



parameter in ([5j) for each choice of the chiral scale fi. The scale invariant LEG is defined by 

,2 



5 = -1927r^/g(^)-log^. 

A* 



and the corresponding convention for the S parameter [3] is 



S 



127r 



1927r2z^(^)+log(^Vm^) 



(7) 



(8) 



III. CONTACT TERMS IN WARD IDENTITIES 



Contact terms in the Ward identities can yield finite contributions to the Fourier transform, 

V/A 

SO that Tlfj^i, is no longer transverse. However the contact terms cancel exactly in the difference 
between the vector-vector and the axial-axial correlator for DWF in the — > oo and massless 
limit, provided the conserved/local correlators are used. Any power-divergent contribution also 
cancels in this difference as a result of the chiral symmetry of DWF. 
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Figure 1: Ward identity violations in the chiral limit for the vector and axial vector currents (left) and their 
difference (right). The same scale is chosen in both plots to visualise the cancellation. 



Following the notation introduced in the previous section we test the vector and axial Ward 

V/A 

identities by extrapolating Qu^fiu to the chiral limit. We find that Ward identity violations 
are very similar for both currents and contact term contributions are greatly suppressed in the 
difference as shown in Figure [ij The cancellation in the chiral limit also hints at only small 

effects from the non-conservation of the axial current due to the residual mass of Domain Wall 
Fermions. The DWF axial current acquires a small multiplicative renormalisation which vanishes 
in the — > cxd limit [HI]. Since our simulation is done at a fixed length in the fifth dimension 
(L<j = 16) we have to consider this effect as part of our error budget. In Refs. |40 | I4 H H2] it has been 
argued that A = \Zj[ — 1| receives contributions from both delocalized modes above the mobility 
edge Ac, and from localized near zero modes of H\y. The former is proportional to 6"^^=^" while 
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Figure 2: Left: cancellation between 11*^^' and XI^'^) as a function of q^, Right: relative error for given A for 
the lowest momenta. 
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the latter is proportional to These correspond to two different contributions to the residual 
mass: the exponential piece is linear in the corresponding component of rrij-es while the latter is 
quadratic in the localised contribution to mres that is larger in our case. 

As a pragmatic approach we vary A between and 3amrcs and estimate the relative error on the 
difference 11^ as A ^^^'^ . Our conclusion is that there is no large cancellation for the local-conserved 
correlators since g^II^^)'^ approaches zero, as shown in Figure [2] We assume a conservative three 
percent systematic error in II'^^) for the non-conservation of the axial current. 

IV. RESULTS 

The data presented is from the ensembles generated by the RBC and UKQCD collaboration with 
the Iwasaki gauge action at /3 = 2.25 which corresponds to a lattice spacing of = 2.33(4) GeV. 
The details of the ensembles will be published in [Uj . We simulate with three values of the light 
quark mass which correspond to pion masses on the range of 290 MeV < 171-^ < 400 MeV. Our 
chiral fits for ^5 rely on previous measurements of the pseudoscalar mass and decay constant at the 
unitary quark masses. For convenience we summarize the results obtained from correlators with 
gauge-fixed wall sources in Table|T]along with the vector and axial vector ground state masses which 
are used as a consistency check using Weinberg's sum rules and our data for l^. In our analysis 



ami 


amTT 




amy 


amA 


0.004 


0.1269(4) 


0.0619(3) 


0.356(6) 


0.522(13) 


0.006 


0.1512(3) 


0.0645(3) 


0.366(5) 


0.543(18) 


0.008 


0.1727(4) 


0.0671(3) 


0.388(6) 


0.551(9) 



Table I: Meson masses and pseudoscalar decay constant in lattice units. 

of n'-^-' we include spatial momenta up to (1,1,0) (or equivalent). We find excellent agreement 
between results with and without spatial momentum for the local-conserved II^^^ Figure [sj This 
is in contrast to the local-local data where there are two distinct branches of points. 

As a first step we obtain effective values for Lio for each mass and momentum from our data for 
n^^), Eq.(|3]). Our results are summarised in Table [ll} Since ChPT is known to have only a small 
radius of convergence in p'^, we restrict ourselves to only the lowest momenta. 

Now we perform a one parameter fit to obtain our central value for Liq. In Figure |4] the 
dependence of Liq on the momentum fit-range is shown for a fit including all three masses. In 
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Figure 3: q^H^^^ at low momentum. The different symbols correspond to momenta of type (n + 1,0,0,0): 
squares, (n, 1, 0, 0): circles, (n, 1, 1, 0): triangles with n = 0, 1, 2 

the lower panel we observe that the X'^/d.o.f. becomes larger than one when more than the lowest 
two momenta are included. Fits with fewer mass values exhibit the same behaviour. We however 
conclude that including only the lowest momentum gives a more reliable result since the higher 
momenta yield values for g^XI^^^ which are consistently below the fit curve when including only 
the smallest momentum. Figure [5] Our central value therefore is obtained from a fit to all three 
masses and the lowest momentum only. In Figure [5] the data points included in the determination 
of Lio and the error band of the resulting fit is shown plotted against mass (left) and momentum 
(right). Our value for L\q is 

L\^{fi = 0.77 GeV) = -0.0057(ll)stat(7)sys. (9) 

The systematic error has several contributions: fit-range, lattice artefacts, scale setting, strange 
quark mass and finite volume, which are described in more detail below. 

We estimate the error on the chiral fit by varying the fit-range in mass and and obtain an error 
of 0.0006, or 11%. The leading lattice artefacts with Domain Wall Fermions are parametrically 
of order o^Aq^jj-,. We assume Aqcd = 300 MeV, use = 2.33(4) GeV and double the result to 
obtain a three percent error. The uncertainty in the lattice scale is relevant for the definition of the 
chiral scale = mp. We vary within its error and obtain an additional two percent error. The 
dynamical strange quark mass in this computation is fixed to am = 0.03 which differs substantially 
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Table II: Results for g^II^^^ and the resulting effective Lio for the lowest momenta and all three masses. 

from the physical value. Using a reweighting procedure for the strange quark determinant |44j we 
vary the strange quark mass down to the physical point in lattice units and find variations in L\q 
of less than three percent. Possible finite volume effects (FVE) are well under control, since at the 
lightest mass we have ninL = 4 and the estimate for FVE for f-^ from [35] is 0.5%. Adding the 
errors in quadrature we find a total systematic error of 0.0007 which is dominated by the error on 
the chiral fit. 

Comparing our result with the previous lattice determination [32] and phenomenological esti- 
mates |46| ITT] we find it to be well consistent, Figure [6j 

Lattice data for the lowest states in the vector and axial vector channels allows us to crosscheck 
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Figure 4: Upper panel: dependence of Liq on the fit-range in (ap)^. Lower panel: reduced of the 
corresponding fits described in the text. The dashed vertical line denotes the chiral scale in lattice units. 
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Figure 5: Data set and fit band used for the final result for L 
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om result with the following sum rules for the spectral densities Pv/A |1H1 US] i 

j dss{pv{s)-pA{s))=0, (10a) 

I ds{pv{s)-pA{s))=fl (10b) 

j ds J {pv{s) - pa{s)) = -8lio. (10c) 
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Figure 6: Comparison of Liq with other determinations. 



We assume that the spectral densities are saturated by the hghtest resonances, 



PV/A- fv/A^i^ -'^V/a)- 



(11) 



The resulting simplified sum rules are solved for Liq using our values for mvjmyi and /^r from 
Table I 



{fvmvf - {fArriAf = 0, 
fv~fA = fw^ 



fy _ /I 

2 2 

rriy m\ 



-8L 



10- 



(12a) 
(12b) 

(12c) 



Thus we obtain an independent estimate for Liq for each quark mass which we convert to L^^^mp) 
using the pions masses from Table |Ij In Figure [7] these results are shown by the black circles and 
are compared with our best estimate ([9|, indicated in the figure by the blue band. We find them to 
be well consistent. This agreement is a non-trivial finding since it is found at finite lattice spacing 



and without a chiral extrapolation for the values obtained from the sum rules (12a 12c). 
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Figure 7: Results for L\q from the sum rules in Eq. (12) These are compared with the result from Eq. (|9| 
which is indicated by the blue band. 



V. PION MASS SPLITTING 

The computation of the vacuum polarisation functions over the complete q'^ range allows the 
computation of the electromagnetic contribution to mass splitting between the charged and neutral 
pions |50j, 

ml^-mlo = Td^z^^ = 1261 MeV2. (13) 

The simplest possible ansatz to describe the mass and momentum dependence of n^^-* which respects 
the sum rules ( 12a|12b ) is 

2-rT(l)<' 2 \ r2 , fv^ 

fv = xi+ X2ml, my = X3 + x^ml, (14) 

JA 

where xi, . . . , X4 are fit parameters. We include all data points up to (ag)^ < 1.0 in the simultaneous 
fit to all three masses and obtain a stable fit with x^/d.o.i. = 1-04. Extrapolating to the chiral limit 
we find 

3a . ,q^nW 



4tt 



/ dg2:i^_ = 1040(220) Me V^. (15) 
Jo / 
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We have also tested the two different fit forms used in [3? which include an additional term, 
which is relevant at high q^, however we do not find a significant change in the fit parameters 
xi, . . . , X4 or a smaller x^/d.o.f. when we include these additional fit parameters. An illustration of 




Figure 8: Left: fit with ansatz (14 1, where the grey band denotes the result in the chiral Hniit. Right: mass 
splitting rn?^i 



^± — TO^o as a function of the cut-off momentum q. 
value, the dashed line marks our chosen cut-off. 



2 

max ' 



the horizontal line denotes the physical 



the dependence on the cut-off is shown in Figure |8j The remaining part of the integral is estimated 
under the assumption that the functional dependence follows the asymptotic behaviour for large 



„2 



n(i) ^ (16) 



Matching (16) to our fit result in the chiral limit we obtain c = —0.0062(58) GeV . Here we have 
doubled the error to allow for possible deviations from the asymptotic form. The remaining integral 
is then 

^2 



3a 
47r 



A dg'^/' = 140(130) MeV2, (17) 
yielding a pion mass splitting of 

m^± - mlo = 1180(260) MeV^. (18) 
This error does not include a rigorous investigation of the systematic uncertainties as was done in 



Sec. IV for Lio, however we expect these uncertainties to be well within the large statistical error 
of twenty percent obtained from the integral (15) with < (ag)^ < 1. The value we obtain is in 



excellent agreement with the experimental value 'rv?_± — m?o = 

1261 MeV^ jSU. 
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Finally we can use the result of the fit defined in Eq. (14) to obtain another independent 



determination for Liq. From the slope of the ansatz (14) in the chiral limit, we find LiQ^nip) = 
5.19(14) • 10^'^ which is in agreement with our ChPT analysis. This is a nice check of the saturation 
of the Weinberg sum rules. Unfortunately the large error on our determination of I5 does not allow 
a more quantitative comparison. 



VI. CONCLUSIONS 



We have computed the S parameter in QCD using the gauge configurations generated by the 
RBC-UKQCD collaboration for 2+1 flavors of dynamical DWF fermions. According to the proce- 
dure outlined in Ref. [32], the S parameter is extracted from the form factor that appears in the 
parametrization of the VV-AA correlator. Our final result is 

L[o{l^ = 0.77 GeV) = -0.0057(ll)stat(7)sys; (19) 

where the error is still dominated by the statistical precision of our simulation. On the other hand, 
the systematic error is dominated by the choice of the fit- range used in the chiral extrapolation. 
A further improvement to these simulations would be to include partially twisted boundary con- 
ditions, enabling access to smaller momenta where ChPT can be employed reliably. Using results 
from meson spectroscopy, we were able to check the saturation of the Weinberg sum rule by the 
lowest-lying resonances in QCD. At the current level of accuracy, we find that the contribution 
from the lowest-lying resonances accounts for the total value of the S parameter, which confirms a 
widely-used assumption in phenomenological studies. Our best estimate for the S parameter ([s]) 
with a Higgs boson mass of mn = 120 GeV is 

5 = 0.42(7), (20) 

where we have rescaled the renormalization scale fi with the ratio of the pion decay constant to 
the Higgs vacuum expectation value v = 246 GeV. 

From the same correlators we have extracted the electromagnetic pion mass splitting, and the 
result Ami = HSO ± 260 MeV^ is in excellent agreement with the experimental results. Despite 
the large error, the present computation shows that it is possible to extract this quantity from the 
currently available DWF ensembles. 

The results presented here are only made possible through the chiral properties of DWF which 
ensure the cancellation of the power-divergent contributions in the current correlators. The QCD 



13 



analysis we have performed lays the groundwork for the computation of the S parameter in potential 
Technicolor model candidates such as those proposed Refs. |13l [52] . Computing the S parameter 
in Technicolor theories is a crucial step in identifying the models that survive the constraints from 
electroweak precision measurements. We plan to extend this computation to the models recently 
simulated in Refs. |24j . 
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